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Abstract. We describe an experimental implementation of a free-space communica- 
tion link encoding information in orbital angular momentum (0AM) spatial modes. 
The free-space path was constructed to represent propagation through an equivalent 
1km horizontal turbulent path with C^-^ = 1.8 x 10^^^ m^^/'^. Turbulence induces a 
spreading in the 0AM spectrum leading to cross-talk in the communication channel. 
Methods of dealing with this loss in channel capacity of the communication channel are 
demonstrated using both adaptive correction of the turbulence as well as optimization 
of the channel encoding. 
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1. Introduction 

It is well known that beams of light that contain a phase vortex of the form ?/'(r, 0) oc 
exp{ii(f)), where i is an integer, carry orbital angular momentum. A single photon 
prepared in such a state will have an angular momentum equal to ih in addition to 
any momentum carried by the polarization [1]. In recent years there has been a good 
deal of excitement based on encoding information onto such 0AM beams of light [2H1] . 
Such spatial mode encoding schemes allow for scaling of the bit rate as each pulse may 
contain more than 2 possible symbols. This scheme also provides an enhancement to 
the security when used for quantum key distribution [5l|6]. 

Such an encoding scheme generally implies communication over a free-space 
channel. A known limitation of free-space communication schemes is degradation of the 
signal due to atmospheric turbulence, and much work has been done to study how this 
turbulence affects communication systems that utilize spatial mode encoding [7HT3]. 
However, much of this work focuses on 'weak' turbulence, where the effect is fully 
described by a random thin phase screen in the receiver aperture. The turbulence obeys 
'Kolmogorov statistics' characterized by a single turbulence strength parameter, D/ro, 
where D is the aperture or beam size and tq is a correlation length scale known as 
Fried's parameter p3|[T5]. 

In more realistic situations, such as communication along a long horizontal path 
through which the turbulence is continuously distributed, one will see amplitude 
fluctuations (or scintillation) that impacts one's ability to communicate in addition 
to the problems caused by pure phase fluctuations. Some of this degradation can 
still be compensated for by phase-only adaptive optics (AO) correcting for low spatial 
order aberrations for a horizontal path of a few kilometers or more [16]. However, 
for moderate scintillation one will begin to see intensity nulls that are associated with 
phase vortices [TTlfTS] . These phase vortices or branch points are known to interfere with 
the performance of AO systems [19], and there is a breakdown in the performance for 
horizontal paths around 5 km due to this effect [20]. For communication systems that 
communicate using 0AM this presents an additional problem as these phase vortices are 
precisely the source of the encoding, and randomly generated vortices introduce errors 
in such a scheme. 

2. Simulating thick turbulence in the laboratory 

To more closely examine the effects that a thick horizontal turbulent path might have on 
OAM-based communication while still allowing information transfer and AO correction, 
we chose to consider a 1km path with aperture sizes of the sender and receiver of 
D = 18.16 cm, at a wavelength of A = 785 nm corresponding to a Fresnel number of 
Nf = D'^/{XL) = 42. In addition, the refractive index structure parameter, which 
quantifies the strength of the refractive index fluctuations due to turbulence at any 
point, was chosen to be = 1.8 x 10^^^ m^^^^ along the path, which represents a 
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typical value for a horizontal path near ground level. 

To allow for such a system to be realized in a laboratory setting one must find 
a way to incorporate the turbulence into the channel, as well as properly scale the 
system down to more manageable length scales. This was done by finding an equivalent 
path containing 2 thin phase screens (section 12. ip that can be represented in the lab 
with spatial light modulators (SLMs). Scaling rules that are invariant under Fresnel 
propagation are detailed in section 12.21 and the experimental setup is discussed in 12.31 

2.1. Two phase screen model 

In order to find an equivalent optical path that accurately represents the horizontal 
turbulent channel, it is necessary that the simulated path faithfully reproduce all 
the relevant statistical properties of the actual channel. The relevant parameters for 
0AM communication are the overall strength of the turbulence, the scintillation or 
amplitude fluctuations, the power fluctuation of the beam due to beam wander and 
aperture clipping, and the random generation of branch points. We can represent these 
4 parameters by using thin Kolmogorov phase screens each with its own value of tq. 
The values of tq for each screen, as well as each screen's position along the path, give 
us four independent parameters that are tuned until the 2 screen path reproduces the 
same 4 parameters of the full horizontal channel. 

The first step in designing a two-phase-screen path to simulate horizontal turbulence 
is to quantify the 4 statistical parameters of a thick path. The overall strength of 
turbulence is represented by D/ro, just as in the thin-phase-turbulence model. Assuming 
a constant value of the refractive index structure parameter along the path, C^z) = C^, 
then this parameter is given by the expression [H] 
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which for our system is found to be Tq = 24.4 mm. A well known measure for scintillation 
is the log-amplitude variance [21], given by the expression 

4 = 0.5Q3k'/' f CI{Z)Z''" dz = 11 ^7/6^2^11/6^ (2) 

which for our system is found to be cr^ = 0.197. The normalized power over the aperture 
S of area As, defined by P = ^ drexp (2x(r)), where is the log-amplitude of 
the field, is used to numerically find the normalized power variance in the aperture: 

where C^(r, r') = ((x(r) — (x)) (xl^O ~ (x))) is the log-amplitude covariance 
function [22] • This value is computed numerically for our system giving ap = 7.04 x 10~^. 
Finally we require the average density of branch points created by the turbulent 
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path, pbp, to be preserved. This density is found by Monte-Carlo simulation to be 
p^p = 500/m2. 

The second step in designing the two-phase-screen path is to find the values for the 
position and tq for each screen that will give the same values of D/vq, a^, ap, and psp of 
the thick path. One starts with an initial guess for roi (i.e. ro for screen one), and then 
solves for the value of ro2 that will give the correct value for D /r^. Next, one randomly 
picks a value for the position of the second screen, 2:2, and finds Zi such that one gets the 
correct value for a^. Z2 is then varied (along with zi to maintain to set the correct 
value for ap. Given this solution, ppp is computed by Monte-Carlo simulation. If at this 
point one gets the correct ppp ^ solution has been found, otherwise one starts over with 
a new choice for tqi. Using this procedure we found we could simulate our 1km path 
with the parameters roi = 3.926 cm, ro2 = 3.503 cm, zi = 171.7 m and Z2 = 1.538 m 
(measured from the sender's aperture). 

1 
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Figure 1. Probability distribution of the Strehl ratio for a thick turbulent path 
represented by 10 Kolmogorov phase screens (red line) and its equivalent 2 phase 
screen solution (blue line). 

As an independent test of this solution, a beam propagation simulation was 
performed to compare the thick turbulent path represented by 10 random Kolmogorov 
phase screens with the analogous 2 screen solution. In each simulation a different random 
realization of turbulence was made and the Strehl ratio, defined as the ratio of the peak 
intensity to ideal peak intensity of a spot at a focal plane of the receiver, was computed. 
By repeating this many times, a probability distribution for the Strehl ratio was found 
and the results are shown in Fig. [1] As can be seen in the plot, the Strehl ratios of the 2 
screen and the 'continuous,' 10 screen paths show very good agreement with each other. 
This demonstrates that the 2 screen model not only reproduces the correct mean values 
for the statistical parameters of interest (by construction), but can also be expected to 
give similar distributions of possible measurement outcomes. 
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2.2. Fresnel scaling 

The second thing that must be done to effectively simulate a turbulent path in a 
laboratory setting is to scale the optical paths down to more manageable lengths. In 
order to ensure that the scaled path still represents the desired physical path, the Fresnel 
equation must remain invariant under the scaling. Fresnel propagation from one plane 



r=(x,y) p=(u,v) 

Figure 2. Propagation from plane A, described by coordinates r — {x, y), to plane B 
with coordinates p — (u, v) 

to another a distance z away as shown in Fig. [2] is given by: 

UBip) = '^ I UA{r)e^^^'^-'^' d'r (4) 

Now if we scale the coordinates using r' = arr, p' = app, and z' = a^z then the above 
propagation equation becomes: 



/ \ / / / \ 2 

f \ t-iTctz I _r p_ 



U^{-]= "\ u I e^V^-S^; dV (5) 

\apj arXz'/a^J \ar J 

Since we require the Fresnel number to remain constant, = ctpa^. Then we can 
rewrite equation [5] as: 

ikz'(l-l/az) / n' \ -np'^ / ^ pikz' r / r' \ i2 . ,, \ 

Ub { ^ ) = e-'^ Ua -] e-'^^e^^'-^'^'dh'] (6) 



\z' J \a. ^ 

where fp = z' / ^1 — and fr = z' / ^1 — . From Eq. E] we see that the horizontal 
path between planes A and B can be scaled down (to within a scaling and phase 
constant) simply by adding a lens with focal length fr at A, propagating a distance 
z', and then adding another lens with focal length fp at B to cancel out the residual 
quadratic phase. 

2.3. Experimental setup 

A diagram of our experimental setup is presented in Fig. [31 The sender, Alice, prepares 
a beam in a specific 0AM state to send to the receiver. Bob, from an attenuated HeNe 
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Figure 3. Alice sends a beam prepared in a specific 0AM state, £, to Bob. Bob 
receives the beam after propagation tlirough a channel representing a 1 km turbulence 
path. The beam is (optionally) corrected using a deformable mirror and sent to a 
sorter to make a measurement of the of the 0AM spectrum of the beam. 



using an SLM combined with a 4/ system as described in [23l|2l]. The prepared state 
is then sent through the simulated 1km path scaled down as described in section [2^ to 
a total length of 1.3 m. The two thin phase screens used to simulate thick turbulence in 
our setup (section 12. ip were implemented using an SLM in a double pass configuration. 
In addition, the quadratic phases required for proper scaling of the propagation path 
were added to the phases on the SLMs. 

After propagation through the turbulent channel, the beam at Bob's aperture is 
imaged with a 4/ system onto a Thorlabs adaptive optics (AO) kit consisting of a 
12 X 12 actuator deformable mirror and a Shack-Hartmann wavefront sensor. After 
the AO system, the beam is similarly imaged onto the first element (Rl) of the 0AM 
sorter. The 0AM sorter uses two refractive optical elements (Rl and R2) and a Fourier 
transforming lens (FT Lens) to spatially separate different 0AM modes allowing for the 
0AM spectrum to be efficiently measured as described in |25] . 

3. Experimental results 

In order to examine the effects of the turbulent channel on 0AM communication, we 
experimentally measured the 0AM spectrum that Bob detects conditioned on what 
Alice sent. In a perfect channel, if Alice sends 0AM mode s, then Bob will measure an 
0AM spectrum that is simply a Kronecker delta centered at the same mode. However, 
in an imperfect or turbulent channel, there will be some spreading into neighboring 
0AM modes to the prepared state. The conditional probability matrix, P{d\s) where 
d is the detected 0AM mode and s is the sent mode, provides a natural expression for 
this crosstalk induced by the imperfections or turbulence in the channel. 

P{d\s) is plotted for 3 different scenarios in Fig. |H Fig. shows P{d\s) when 
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a) No Turbulence b) Turbulence With No Correction c) Turbulence With AO Correction 




Detected Mode Detected Mode Detected Mode 

Figure 4. Measurement of the crosstalk in tlie cliannel represented by tire conditional 
probability matrix, P(d\s) for three cases: a) no turbulence, b) with turbulence, and 
c) with turbulence and adaptive correction. 

there is no turbulence in the channel, showing only crosstalk due to any misalignment 
in the system and inherent crosstalk of the sorter [25]. Fig. |3)d shows the effects of 
thick turbulence ensemble averaged over 100 realizations, which act to greatly spread 
the signal over many neighboring channels. This selective spreading of 0AM into 
neighboring modes rather than randomly into any 0AM state is qualitatively similar 
to what is seen in the thin turbulence regime demonstrated in [9]. For the third case 
shown in Fig. Ht, adaptive correction was applied to the turbulence with the AO system, 
allowing much of the signal to be recovered. The phase aberrations induced from each 
realization of turbulence was sensed and corrected by the AO using the 0AM i = 
mode. Each mode was then sent through the channel and AO system and the 0AM 
spectrum was measured by Bob. This procedure was repeated and averaged over 50 
realizations of turbulence. 
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Figure 5. Channel capacity as a function of the number of 0AM modes used for the 
spatial encoding. 

In order to quantify the crosstalk induced by the turbulence as well as the quality 
of the AO correction, we compute the mutual information between Alice and Bob for all 
three cases above. The mutual information between Alice and Bob, I{A;B), provides 



Optical 0AM communication through 1 km horizontal path 



8 



a measure of the channel capacity, as it gives the maximum possible transmission rate 
that can be extracted by Bob in bits per photon (or classically per detected pulse). The 
mutual information is given by the expression, 

where is number of distinct symbols. Fig. Hlqualitatively showed that thick turbulence 
greatly impacts the quality of the channel. Using Equation [TJ we quantify these results 
by calculating the mutual information as a function of the encoding dimension N. The 
mutual information for the three cases of no turbulence, thick turbulence, and turbulence 
with AO correction, is plotted in Fig. [5] as a function of A^. One can see that the AO 
system allows us to cancel roughly half of the loss of channel capacity due to turbulence. 
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Figure 6. Channel capacity as a function of the spacing between 0AM modes used 
for communication. 

Further, since turbulence preferentially scatters power into neighboring 0AM modes 
rather than randomly into all modes, one can increase the channel capacity by choosing 
to use a less dense set of 0AM modes [11]. For instance, rather than encoding every 
0AM mode, one can encode in every other mode (or third or fourth mode etc.). 
Changing the encoding is also independent of any AO system one may use, and thus a 
modified encoding can be used along with AO correction to further enhance the channel 
capacity. Fig. |6] shows the increase in the mutual information one can obtain for a 
given number of encoded modes, N. The channel capacity of an ideal 2-bit system is 
shown for reference. It is worth noting that the use of spatial mode encoding shows 
an improvement over such a system with very moderate resources (i.e. 3 modes with a 
channel spacing > 4). 

4. Conclusions 

In this work we have demonstrated how one could experimentally simulate a thick 
horizontal turbulent path as part of a free-space communication channel. This model 
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was applied to the case of a 1 km path and the effects were studied in the context of 
an 0AM communications channel. It was shown that although the turbulence severely 
degraded the channel, a high information capacity was still possible by optimizing the 
encoding and adaptively correcting for some of the induced phase aberrations. 

We acknowledge Miles Padgett, Martin Lavery, and Daniel Gauthier for helpful 
discussions. Our work was supported by the Defense Advanced Research Projects 
Agency (DARPA) InPho program. 

References 

[1] Allen L, Beijersbergen M, Spreeuw R and Woerdman J P 1992 Phys. Rev. A 45 8185-8189 

[2] Gibson G, Courtial J, Padgett M J, Vasnetsov M, Pas'ko V, Barnett S M and Franke-Arnold S 

2004 Opt. Express 12 5448 
[3] Wang J, Yang J Y, Fazal I M, Ahmed N, Yan Y, Huang H, Ren Y, Yue Y, Dolinar S, Tur M and 

Wilner A E 2012 Nat. Photonics 6 488-496 
[4] Boyd R W, Jha A K, Malik M, O'Sullivan C, Rodenburg B and Gauthier D J 2011 Proceedings of 

SPIE, vol. 7948 pp. 79480L-79480L-6 (SPIE) 
[5] Bourennane M, Karlsson A, Bjork G, Gisin N and Cerf N 2002 J. Phys. A: Math. Gen. 35 

10065-10076 

[6] Cerf N, Bourennane M, Karlsson A and Gisin N 2002 Phys. Rev. Lett. 88 

[7] Paterson C 2005 Phys. Rev. Lett. 94 

[8] Tyler G A and Boyd R W 2009 Opt. Lett. 34 142 

[9] Rodenburg B, Lavery M P J, Malik M, OSullivan M N, Mirhosseini M, Robertson D J, Padgett 
M J and Boyd R W 2012 Opt. Lett. 37 3735 
[10] Gbur G and Tyson R K 2008 J. Opt. Soc. Am. A 25 225 

[11] Malik M, OSullivan M N, Rodenburg B, Mirhosseini M, Leach J, Lavery M P J, Padgett M J and 

Boyd R W 2012 Opt. Express 20 13195 11204. 57811 
[12] Roux F 2011 Phys. Rev. A 83 

[13] Boyd R W, Rodenburg B, Mirhosseini M and Barnett S M 2011 Opt. Express 19 18310 
[14] Fried D L 1965 J. Opt. Soc. Am. 55 1427 
[15] Fried D L 1968 J. Opt. Soc. Am. 58 961 

[16] Levine B M, Martinsen E A, Wirth A, Jankevics A, Toledo-Quinones M, Landers F and Bruno 

T L 1998 Applied Opt. 37 4553 
[17] Fried D L and Vaughn J L 1992 Applied Opt. 31 2865 
[18] Tyler G A 2000 J. Opt. Soc. Am. A 17 1828 
[19] Fried D L 1998 J. Opt. Soc. Am. A 15 2759 

[20] Primmerman C A, Price T R, Humphreys R A, ZoUars B G, Barclay H T and Herrmann J 1995 
Applied Opt. 34 2081 

[21] Tatarski V I 1961 Wave Propogation in a Turbulent Medium (McGraw-Hill Books) 

[22] Andrews L C and Phillips R L 2005 Laser Beam Propagation Through Random Media 2nd ed. 

Press Monographs (SPIE Press) 
[23] Arrizon V, Ruiz U, Carrada R and Gonzalez L A 2007 J. Opt. Soc. Am. A 24 3500 
[24] Gruneisen M T, Miller W A, Dymale R G and Sweiti A M 2008 Applied Opt. 47 A32 
[25] Lavery M P J, Robertson D J, Berkhout G, Love G D, Padgett M J and Courtial J 2012 Opt. 

Express 20 2110 



